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TECHNICAL PAPER 


DEVELOPMENT OF HOMOTOPY ALGORITHMS FOR FIXED-ORDER MIXED 

H 2 /H m CONTROLLER SYNTHESIS 

I. INTRODUCTION 


Modern control theory has revolutionized control system design over the past decade. H 2 and 
H„ methods have gained widespread recognition and are used in controller synthesis for SISO and 
MIMO problems. A significant disadvantage of modern control techniques is that the resulting com- 
pensator is the same order as the generalized plant, which is often larger than the original plant due 
to the inclusion of frequency-dependent weights to achieve the desired performance and robustness 
characteristics. Real-time implementation of these controllers for large-order systems is prohibitive 
due to the computational burden of fast throughput times. 

One indirect approach to alleviating this is to reduce the order of the plant and to synthesize a 
controller based on this reduced-order design plant. An alternate indirect approach is to design a 
full-order controller and then to apply model reduction to the controller. In either case, indirect 
methods cannot guarantee closed-loop stability and are suboptimal in performance. Direct methods 
impose constraints on controller order or architecture. In an optimization based synthesis procedure, 
necessary conditions are formulated for the constrained closed-loop system that ensure internal 
stability. The optimal projection approach of reference 1 is an example whereby order constraints are 
imposed on the controller, and the necessary conditions for minimizing a quadratic cost functional 
with respect to the fixed-order controller are derived. The resulting necessary conditions consist of 
two modified Riccati equations and two modified Lyapunov equations coupled by an oblique pro- 
jection matrix. However, solution of the necessary conditions for realistic large order systems is a 
difficult task. Reference 2 employs homotopy methods to solve the optimal projection equations. 
Newton methods have also been applied. 3 

Optimal projection has also been extended to LQG control with an //„ norm over-bound. 4 This 
mixed // 2///00 optimization problem seeks to minimize the H 2 norm of one transfer function for perfor- 
mance while satisfying a bound on the //<*> norm of another transfer function for robustness. The true 
mixed problem has two inputs and two outputs, indicating different classes of disturbances and per- 
formance variables. Much of the research into the mixed problem considers variations of the true 
problem with only one input or only one output. Reference 4 considered the case ot two outputs and 
one input, with both full-order and fixed-order control. The difficulty here is the size of the gap 
between the H M over-bound and the true H norm. 5 Reference 6 addresses the two input, two output 
problem with output feedback including the fixed-order problem, but does not attempt to solve it. 
Recently, reference 7 used a differential game formulation to obtain fixed-order controllers for the 
true mixed problem. A conjugate gradient technique was applied to solve these resulting necessary 
conditions. 

The objective of this paper is to build on the results of references 4 and 7 by presenting 
homotopy algorithms for solving the H 2 , and true mixed // 2 ///„ fixed-order compensator syn- 
thesis problems. The paper is organized as follows. Section II presents a formulation ot the problem 
with the compensator in controller canonical form. The necessary conditions for the fixed-order H 2 
controller are developed. These results are then extended to the fixed-order H„ and mixed H 2 tH^ 



compensator design using the differential game results of reference 7. Section III introduces homo- 
topy methods and develops homotopy algorithms for solution of the H 2 , H„, and the mixed 
problem formulations of section II. Section IV presents numerical evaluations of these homotopy 
algorithms, and section V concludes the paper. 


II. PROBLEM FORMULATION 


The generalized plant of a standard control problem is given by: 

x = Ax+BiW+B 2 u , ( 1 ) 

z = C l x^D n u , (2) 

y = C2X+D21W+D22U , ( 3 ) 

where x e R n is the state vector, we R nw is the disturbance vector, u e R nu is the control vector, 
z e R nz is the performance vector, and y e R ny is the observation vector. It is assumed that: 

• (A.Bj.Cj) is stabilizable and detectable 

• CA.B2.C2) is stabilizable and detectable 

• D n has full column rank 

• £> 21 has full row rank. 

A general compensator for this system is 

x c — A c x c +B c y , (4) 

u = C c x c , (5) 

where x e R nc is the state vector of the controller, the dimension of which can be specified. Figure 1 
illustrates this design framework. Closing the loop using negative feedback yields the closed-loop 
system dynamics: 


* = Ax+Bw , 

z = Cj , 


where 


x = 



(6) 

(7) 


( 8 ) 


2 



(9) 


A=\ A 

B c C 2 


-B 2 C c 

ArB<PtfCc 


B = 


Bi 

BJ>2 1 ’ 


( 10 ) 


r - [Cj -D l2 C c ] . 


( 11 ) 



Figure 1 . Generalized plant with general compensator. 

The set of all internally stabilizing compensators is defined as: 

S c = { ( A c ,B c ,C c ) : A is asymptotically stable } . (12) 

For an H 2 problem, the objective is to minimize the // 2 -norm on the closed-loop transfer 
function from disturbance inputs to performance outputs: 

T zw = €(sI-A) l B , (13) 

where the disturbances are confined to the set of signals with bounded power and fixed spectra. This 
leads to three equivalent H 2 optimization problems. For impulsive inputs Wj = 8{t), the objective is 

imn|j(^,S c ,C c ) = || r w || 2 = (S||z,||) 1/! , i=l ny ) , (14) 

where zi is the response to w, and 1 1 z- t \ | 2 denotes the L 2 norm. For w,- = ai5(t) with £{a,-}= 0, 
E{aidj) = 5(i-j) the objective becomes: 

min ^J(A c ,B c ,C c ) = E M z(t) T z(t)dt^j , (15) 

where z(t) is the response to an initial condition x(0), and E x0 { • } denotes the expectation over a 
distribution of initial conditions defined by E{x( 0)} = 0, E{jc(0)jc(0) r } = SB T . If the disturbance is 
modeled as white noise, the objective is: 
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(16) 


min |7(4,B c ,C c ) = ^lim^ E{z(t) T z(t)}} . 

It can be shown that the cost in the three formulations given above can be expressed as: 

J{A c £ c ,C c ) = tr{QB$ r ) = tr{PC T C) , (17) 

where 

AP+PA t +SS t = 0 , (18) 

A T Q+QA+C T C = 0 . (19) 


P is the controllability grammian of (A,S) and Q is the observability grammian of (C, A). 

In order to obtain the # 2 -optimal compensator, the Lagrangian is defined as: 

£(Q,LA c ft c ,C c ) = tr{QBB T HA T Q+QA+C T C)L} , (20) 

where L is a symmetric matrix of multipliers. Matrix gradients are taken to determine the first-order 
necessary conditions: 


d£ djC d£ d£ djC 

dQ = dL"dA c = W c "dC c = {) ' (21) 

Hence, computing an ^-optimal controller of fixed-order nc < n for the general controller structure 
given in equations (4) and (5) requires the simultaneous solution of five coupled equations. This is 
not only computationally expensive, but is also further complicated by the fact that the problem is 
overparametrized with such a compensator. 

To avoid the problem of overparametrization, a canonical form description for the controller 
can be used. It was shown in reference 8 that if either a controller or observer canonical form is 
imposed on the compensator structure, the number of parameters is reduced to its minimal number. 
The internal structure of the compensator is prespecified by assigning a set of feedback invariant 
indices v,. In controller canonical form the compensator is defined as: 


x c - P°x c +!\f ) u c -t/ > y , (22) 

u c = -Px c , (23) 

u = -Hx c , (24) 

where x c e R nc and u c e R ny . P and H are free-parameter matrices, and P° and N° are fixed matrices 
of zeros and ones determined by the choice of controllability indices v, as follows: 

P° = block diag^ 0 ,...,/*} , (25) 
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N° = block diag{[0 ... 01]^ Xv J , 

where /= 1, ny. The controllability indices must satisfy the following condition: 

2j v / = nc, i=l, ...,ny . 

i 


( 26 ) 

(27) 

(28) 


Figure 2 shows the structure of the controller. Similarly, a compensator in observer canonical form 
can be constructed. In this paper, only the controller canonical form is employed, which imposes the 
lower bound nc > ny on the order of the compensator. 



(29) 




(30) 

(31) 

(32) 

(33) 


Equations (30) through (33) define a static gain output feedback problem where the compensator is 
represented by a minimal number of free parameters in the design matrix, G. The augmented system 
is shown in figure 3. The closed-loop system is given by: 
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* = ( A-B 2 GC 2 )X+B x w = Ax+Bw , 
z = (C\-D X1 GC 1 )X = Cx . 


( 34 ) 



(35) 


Figure 3. Augmented system with compensator in controller canonical form. 

Minimizing the // 2 -norm of T zw = C(sI-A)~ l B utilizes the same Lagrangian as given in 
equation (20), but now £ is only a function of three parameter matrices, i.e., £{Q,L,G). Thus, only 
three first-order necessary conditions result: 


30 = Al+lA t +BB t = 0 , (36) 

gf = A T C+a?+C r C = 0 , (37) 


OJu t T T T 

— = 2(D n D n GC 2 -D n C l -B 2 Q)LC 2 =0 . (38) 

3G 

This demonstrates that using the controller canonical form yields simpler expressions for the neces- 
sary conditions, with the additional benefit of minimizing the number of compensator parameters. 

Controller canonical forms can also be used to solve the H problem. The objective is now to 
minimize the 00 -norm of the transfer function from disturbance inputs w to performance outputs z 
given in equation (13). In this case, the necessary conditions for an //« -optimal fixed-order com- 
pensator gain G are: 7 


!£ = a + r 2 BB T QjL + La+r 1 BB T Qj T + BB T = 0 . o9> 

g£ = & 1 Q.,+Q.J+e T £+Y- 1 Q.M T Q..,=<> . (40) 

— = l{D\ 2 D n GC 2 -D\ 2 C x -B 2 QJLC T 2 =Q , (41) 
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where 

£«2»X,G) = tr{ Q„M T + (A T Q^+Q^A+C T C+ y 2 QjffQJL } . (42) 

As in the H 2 problem, three coupled equations have to be solved to obtain a fixed-order compensator 
which satisfies the constraint HT^H.. < y. 

Using this approach, fixed-order 7/oo-design can also be extended to fixed-order ^/-synthesis. 
Since H „ controller design is a subproblem when designing for robust performance with structured 
uncertainty, the fixed-order technique introduced above has the potential to constrain the order of the 
controller which is normally subject to significant increases in the /r-synthesis procedure. 

The mixed H 2 /H ao problem can be approached in a similar fashion. In this case, the gener- 
alized plant has additional inputs and outputs w p and z p , respectively, which define the H 2 perfor- 
mance criterion (fig. 4). The inputs w and outputs z are used to define the //« performance criterion. 
Using the controller canonical form for the compensator, the augmented system for the mixed prob- 
lem is: 


where 


i = Ax+B\\v+B p w p +B2U , 
z p — CpX+DipU , 
z = C]T+D > 

y = c 2 x , 

17 = -Gy , 


B P = 


B P 

-N°D, 


2pJ 

Cp = [Cp 0] , 
D lp =[D lp 0] . 


A 

B 

p 

B 

i 

B 

2 

C 

p 

0 

0 

D 

ip 

c 

1 

0 

0 

D 

12 

c 

D 

D 

D 

2 

?£_ 

21 

22 


Figure 4. Mixed H 2 !H problem. 


(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 
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The other expressions are the same as in equations (30) through (33). Consequently, the closed- 
loop system is given by: 


x = ( A-B 2 GC 2 )x+B p w p -\-B x w = Aj+B p w p +Bw , 

(51) 

z p = (C p -D Xp GC 2 )x = C p * , 

(52) 

z = (C X -D X2 GC 2 )J = Cj . 

(53) 

In order to formulate the performance index of the mixed problem, the Lagrangian for the H 2 problem 
in equation (20) is adjoined to the Lagrangian for the H «, problem in equation (42) by a scalar weight 
A: 

£ = tr{ Q m BB T HA T Q^+Q^A+C T C+y~ 2 Q x BB T Q M )L+Ax€ T p C p +(AX+XA T +B p B T p )L p } . (54) 

The weight A on the // 2 -norm allows a tradeoff between (H 2 ) performance and the H„ 
order necessary conditions are: 

norm. The first 

aOo = ^ + r -2 ^ r a»)L+L(A+y- 2 5s r a 0 ) r +fiB 7 '= 0 , 

(55) 

§f = Z T Q„*<>.,A+e T £+y- 1 Q,MQ. = 0 , 

(56) 


(57) 

~Ax + xA\B p B T p = 0 . 

(58) 

^ = 2(d] 2 d X2 gc 2 lcI-d] 2 c x lcI-bIqj.cI+adI p d 1p gc 2 xcI 


-Ad\ p c p xc t 2 -b\l p XC T 2 ) = 0 . 

(59) 


As demonstrated above, imposing a controller canonical form on the compensator structure 
provides a powerful tool for the design of fixed-order controllers. Promising results have been 
obtained for the and the mixed problem in reference 7 where a conjugate gradient method was 
used in the computation. A disadvantage of this method is that convergence slows down near the 
optimum. Also, an initial starting guess for the compensator gain G has to be provided that stabilizes 
the closed-loop system. In this paper, a homotopy method is used to continuously deform the 
solution of a simple problem formulation to the solution of the desired problem formulation. 
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III. HOMOTOPY METHODS 


Homotopy methods offer an attractive alternative to more standard approaches of optimal 
controller synthesis such as sequential and conjugate gradient methods. The basic philosophy of 
homotopy methods is to deform a problem which is relatively easily solved into the problem for which 
a solution is desired. 

Homotopy (or continuation) methods, arising from algebraic and differential topology, embed 
a given problem in a parameterized family of problems. More specifically, consider sets U and Y e SR" 
and a mapping F : U — » Y, where solutions of the problem 

F{u) = 0 (60) 

are desired with ue U and F(u)e Y. The homotopy function is defined by the mapping 
H:Ux[0, 1] 9?" such that: 


H(u l ,l) = F(u) , (61) 

and there exists a known (or easily calculated) solution, u 0 , such that: 

H(u 0 , 0) = 0 . (62) 

The homotopy function is a continuously differentiable function given by: 

H(u{a),a) = 0, Vote [0,1] . (63) 

Thus, the homotopy begins with a simple problem with a known solution (equation (62)) which is 
deformed by continuously varying the parameter until the solution of the original problem (equation 
(60)) is obtained. 10 The power of homotopy methods is that minimization is not strongly dependent 
on starting solution, but depends on local, small variations in the solution. Theoretically, these 
methods are globally convergent for a wide range of complex optimization problems, but in actuality, 
finite wordlength computation often introduces numerical ill-conditioning resulting in difficulties with 
convergence. In light of these numerical limitations, a judicious choice of the initial problem is neces- 
sary for convergence and efficient computation. However, the ability to select an initial problem with 
a simple solution renders homotopy methods more widely applicable than sequential- or gradient- 
based methods, which have a stringent requirement for an initial stabilizing solution. 

Both discrete and continuous methods are used to solve the homotopy. Discrete methods 
simply partition the interval [0,1] to obtain a finite chain of problems: 

H(u,a n ) = 0, 0 = a () < a, < . . . < a N = 1 . (64) 

Starling with a known solution at a„, the solution for H(u,a n+ \ ) is computed by a local iteration 
scheme. Continuous methods involve integration of Davidenko’s differential equation, which is 
obtained by differentiating equation (63) with respect to a, yielding: 
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-1 


dH 
da ' 


( 65 ) 


du _ 
da~ 



Given w(0) = «o, this initial value problem may be numerically integrated to obtain the solution at 
a = 1 if the solution exists and is uniquely defined. 

Research remains to be done in the application of homotopy algorithms. Efficient application 
of homotopy methods depends on the initial problem, the final problem, and the deformation under- 
taken. Given a good initial solution, the key to convergence is the ability to accurately track the 
solution curve, which is determined by the deformation undertaken. The ability to predict the solution 
along the homotopy path via Davidenko’s differential equation makes continuous homotopy methods 
superior to discrete methods. Efficient computation of the Hessian is the primary issue for practical 
implementation of continuous homotopy algorithms. In the following sections, continuous homotopy 
algorithms are presented for reduced-order H 2 ,H„, and mixed H^H^ compensator designs. 


A. Homotopy Algorithm 

This section describes the algorithm used for implementing the continuous homotopy. In 
essence, a mixed discrete and continuous approach is employed where Davidenko’s differential 
equation (65) is integrated along the homotopy path, and at discrete points along the trajectory, a 
Newton’s correction is used for local optimization to remove integration error. Newton’s method, 
which has quadratic convergence properties in a neighborhood of the local minimum, allows a crude 
integration procedure with large step sizes to be employed for efficiently tracking the solution curve. 
This approach follows closely that of references 1 1 and 12 and is employed in the following algorithm. 

1. Find initial solution (a = 0). 

2. Advance a. 


«!,* = «0 + AOo,* . 


3. Predict 6. 


6(a l k ) = e(a 0 ) + Aa 0k O'(a 0 ) , 


where 


*<«>-£- 



dH 
da ’ 


4. Check prediction error. 


e k (6,a) = ||./ e (0(a 1 ,*))|| < e . 
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a. If error less than tolerance, continue. 

b. If not, 0.5Aao k — > Aocq * +1 . 



c. Increment k and repeat steps 2 to 4. 

5. Correct with Newton’s method to compute local minimum. 

6. If a = 1, stop. Otherwise, go to step 2. 

Various approaches may be taken when selecting the deformation, but the general procedure 
applied in this effort is outlined as follows: 

• Synthesize a low-authority H 2 (full-order) compensator 

• Reduce the compensator to desired order and transform to canonical form 9 

• Set /large enough so that the H 2 and compensators are approximately equivalent 

• Use homotopy to deform the initial low-authority, reduced-order H 2 compensator into a 
full-authority reduced-order H 2 compensator (H 2 homotopy) 

• Deform the full-authority H 2 compensator into a nearly optimal //„ compensator with / 
approaching its infimum (H<*> homotopy) 

• At discrete values of A, fix /and deform the compensator into the mixed H 2 IH^ compen- 
sator by varying A (mixed H 2 IH<*> homotopy). 

This procedure was chosen because it has been observed numerically that order reduction 
techniques tend to work best for low-authority LQG controllers. 11 A canonical form is imposed on 
the compensator structure to minimize the number of free parameters, which in some cases can also 
lead to numerical ill-conditioning. A balancing transformation which does not affect the controller 
characteristics relaxes the strict structure in the P° and N° matrices in equations (25) through (27) 
and improves the conditioning of the problem. 

The procedure outlined above separates the compensator synthesis into distinct phases. The 
initial reduced-order, full-authority compensator is synthesized using the H 2 homotopy, which is 
then deformed into the reduced-order H ^ compensator. During the H phase, the scalar H 2 norm 
weight A is fixed (as are the plant matrices) and only the norm overbound / is varied. At 
discrete values, /is fixed and A is varied to perform the H 2 norm minimization. Thus, the procedure 
alternates between the //„ and H 2 norm minimization. 

During the homotopy, both the predicted and corrected gains are checked to assure closed- 
loop stability. After each correction step, the cost gradient is checked to verify descent. During the 
H„, homotopy, the solvability of the Riccati equation using predicted or corrected gains must also be 
checked. If any of these conditions are violated during correction, the correction step size is scaled 
and the condition is checked again. If scaling the correction step size is ineffective, the prediction 
step size is decreased and the prediction phase is repeated. This process continues until the homo- 
topy is completed or until the prediction step size is decreased below a prespecified tolerance. The 
following sections detail the derivations employed in the homotopy algorithm. 
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B. The H 2 Case 


In reference 11, a continuous homotopy algorithm is presented for H 2 compensator synthesis 
when the compensator has a general architecture, which requires the solution of the five coupled 
equations given by equation (21). The following development parallels that of reference 1 1 except in 
this formulation, a controller canonical form is employed for the compensator dynamics which results 
in the three necessary conditions (equations (36) through (38)). 

Consider equations (34) and (35). The necessary conditions for an H 2 optimal compensator 
are given by equations (36) through (38). Define 0 to be a vector comprised of the free compensator 
parameters 


0=vec(G) , (66) 

where G is the output feedback gain matrix defined in equation (33). The gradient of the cost is: 


M = de = vec 


(d£\ n 

(acj - 0 ’ 


(67) 


where d£/BG is given by equation (38). 

The homotopy function is defined as 


d£(0,a) (d£(6,cc)) „ 

me, a) = M = «c(— jjj J=o 


( 68 ) 


Note that £ is now a function of the homotopy parameter a since the system matrices are now func- 
tions of a. The gradient of the homotopy function is 


V[//'(0,a)] = [//„//„] . 

1. Computation of Hessian. The hessian of the cost, Hg, is given by: 

fdH(0,ar 


(69) 


H e = vec 


V 


ma)J . 

H0-)’ J = 




(70) 


where N is the number of free parameters. Using equation (38), 

07/ (0,C(') 0 T T T T 

[KD^D^C^D^Ci-B^LCi] , 


00 , 00 , 

= 2(Dl 2 D l2 G {J) C 2 -BlQ (J> )Lcl+2(Dl 2 D l 2GC 2 -Dl 2 C l -BlQ)L^cl 


(71) 

(72) 


where 


(*) 


u)_m 

- 00, • 


( 73 ) 
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For 6j =g ik . 


G 


u) 

"90, 




lit 


(74) 


which is a matrix of zeros except for a one in the ik element. From equations (36) and (37), deriva- 
tives of L and Q are obtained by solving the Lyapunov equations: 


AL { j \ L { j ) A T +[ A { J \+ LA < J ) T + m T ) {i) ] = 0 , (75) 

A T c! J) +d J) A+[^ j)T QH 2A (J) HC T C) iJ} ] = 0 , (76) 


where from equations (34) and (35), 

A (S) = - B 2 d J ) c 2 , 

(M r ) w = o , 

(€ T C ) iJ) = (C 2 V ) 7 'd[ 2 d 12 gc 2 )+(c[g°' ) 7 'd[ 2 d 12 gc 2 ) 7 ' 

~(c\ D n (f ]) C 2 )-{C] D n C V C 2 ) T . 


(77) 

(78) 


(79) 


2. Computation ofH a . The derivative of the homotopy function with respect to the homotopy 
parameter, a, is: 


H a = vec 


(dH(6,a) 
[ da 


where 


dHC6,a) d r — f — 


= — -[2(d[ 2 D 12 GC 2 -D 12 Ci-6 2 G)LC^] , 
da da 

= 2[(/) 12 D 12 GC 2 +/) 12 Dj 2 GC 2 +D 12 /)] 2 GC 2 -/) 12 C 1 -D] 2 C j -B 2 q 


(80) 


(81) 


-/?[(2)LC 2 +(D| 2 D 12 GC 2 -/>| 2 C,-/i 2 0(LC 2 +LC 2 )] , 


— T — T — — — 7 — — T 

1 . i T ^ 1 i rx 1 Z 1 I > J , 




(82) 


and 


(*) = 


dn 

da 


(83) 


The derivative expressions (equations (80) through (82)) depend on the deformation under- 
taken in the specified problem, i.e., the initial and final problem. In general, suppose that the defor- 
mation of the matrix A is prescibed to be: 


13 



A(a) = £ 0 (a)+a(ifyaM 0 (a)) , (84) 

where the 0 and / subscripts indicate the initial and final system matrices, respectively. It follows 
that 


A=A f -A 0 . (85) 


The derivative of other plant matrices are determined accordingly. The derivatives of L and Q with 
respect to a are obtained from equations (36) and (37) by solving the Lyapunov equations: 


0 = AL +La t + (Al+lA t +§B t +b^ t ) , 

(86) 

0 = A t &QAhA t &QA+C t C+C t C) , 

(87) 

where from equations (34) and (35) 


A — A—B 2 G C 2~B 2 (j C 2 » 

(88) 


(89) 

t = t x -t l2 GC 2 -D l2 Gt 2 . 

(90) 


By employing canonical compensator formulations of equations (30) through (33), the 
expressions for the derivatives of the augmented system matrices reduce to: 


A = 


Af -A 0 0 

-N°(C Zf -C 2fi ) Oj ’ 


(91) 


B l,f~B x0 

-N°(D 2l f -D 2 x 0 ) • 


(92) 


B,= 


^ 2 ,/ - 5 2,0 0 
-N°(D 2 2,r D n,d> 0 


C\ - [ Cy j-Ci Q 0 ] , 

t 2 = [ 0 0] , 

D \1 - [ ^ 12 / - ^ 12,0 0 ] • 


(93) 

(94) 

(95) 

(96) 


Thus, the use of canonical forms not only simplifies the necessary conditions by grouping all the free 
compensator parameters into one feedback gain matrix, but also simplifies the derivative expres- 
sions. The presence of the zero subblocks significantly enhances the computational efficiency of this 
approach. 
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When implementing the procedure described at the beginning of this section, the above 
equations may be further specialized. In this procedure, the initial and final plant matrices are the 
same, and the homotopy is performed only on the measurement and process noise intensities, £>12 

and £> 21 - Hence, A , B 2 ,C X , and C 2 are identically zero. 


C. The Hoo Case 


A continuous homotopy algorithm has also been developed to solve the fixed-order H con- 
trol problem. The development for the //«, homotopy algorithm is identical to that of the previous sec- 
tion for the H 2 homotopy algorithm with the exception of additional terms in the necessary conditions 
(equations (39) through (41)) resulting from the / weighted disturbance term in the cost function 
(equation (42)). Only the portions that differ from the previous section will be presented here. 

As with the H 2 case, when implementing the procedure described at the beginning of this 
section, the previous equations may be further specialized for the //<*> case. In this procedure, the 
homotopy begins with a full-authority H 2 compensator (obtained by the H 2 homotopy) which is then 
deformed into a full-authority compensator by decreasing /. Note that for large values of /, the 
H ^ necessary conditions are equivalent to the H 2 case. The initial and final plant matrices and 
weights are identical, and the homotopy is performed on / only, resulting in considerable simplifi- 
cations in the computation of H a . The value of /is linearly varied from an initial high value toward a 
lower bound (determined from a full-order design) according to: 

Y~ ymax-^/niax-yniin)- (97) 

Thus, the homotopy function defined by equations (68) and (41) is implicitly a function of a. 

1. Computation of Hessian. The hessian for the //„ homotopy is identical to the hessian for 
the H 2 case, except that the observability grammian, Q, is replaced by Qoo, which is the solution of 
equation (39). The hessian is given by: 


^^ = -f{ 2 <D] 2 D l2 GC 2 -D T l2 C r Bl&JLcl\ . 


30, 


33 


= 2(D T l2 D l2 G (J) C 2 -B T 2 Q!d>)L C[+ 2 (D [ 2 D 12 G C 2 -D \ 2 C x -B T 2 QJL^cl 


To obtain expressions for £ (y) and Q^P, differentiate equations (39) and (40) to obtain: 

0 = (A+y~ 2 BB t Q m )L U) +L iJ \A+y- 2 M T Q„) T H(A^+Y~ 2 BB T Q^)L 
+L(A {j) +Y~ 2 BB T Q ^) T ] , 

0 = + QutQjPH&CP) , 


(98) 

(99) 


(100) 

( 101 ) 


where A (}) and (C T t) {J) are given by equations (77) and (79), respectively. 
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2. Computation of H a . The derivative of the homotopy function with respect to the homotopy 
parameter, a, is given by equations (80) and (81) where Q is replaced by Q„, and equation (82) 
reduces to: 


— T — — — T — _ T - _ T _ T t _T 

— = 2[(Z> n G C 2 ~D [ 2 C x -B 2 QJL Cl-Bl&j, c[] . (102) 

L and are obtained by differentiating equations (39) and (40) and are given by: 

o = (A+y^B&q^L +l (A+y~ 2 BB r Q m ) T +(rL+Lr T ) , (103) 

0 = (A+y- 2 BB r Q^ T Q„+<^(A+y- 2 BB T QJA2y- 3 yQjB T Qj , ( 104 ) 

where 

T = -2/“ 3 yBB r Q eo +y- 2 BB T (L . (105) 

The derivative of /with respect to a is obtained from equation (97): 

7 — (/max - )min) • (106) 

Note that a general homotopic deformation of the plant matrices and weights is allowable by intro- 
ducing the variation of the plant matrices with respect to the homotopy. The two-step procedure 
presented here can, in principle, be reduced to a one-step homotopy while simultaneously deforming 
the system matrices (the H 2 case) and /(the case). However, the minimum achievable value of 
/ when designing a reduced-order compensator is usually not known a priori, whereas the 
system matrices must be fully deformed (i.e., a must attain the value of 1) to obtain the desired 
system representation. Thus, the two-step procedure described previously was used in computing 
the fixed-order H „ compensators. For the H„ homotopy, only the overbound / is deformed, thereby 
greatly simplifying the required computation. 


D. The Mixed H 2 IHoo Case 

1. Computation of Hessian. A homotopy algorithm that solves the necessary conditions for 
the mixed H 2 !H aa case equations (55) through (59) is obtained in a straightforward extension of the 
H 2 and Hoc homotopy algorithms. The hessian for the mixed case is given by: 

[^i2^i2GC 2 -Dl 2 C r BlQJLcl^ D T lp D lp GC 2 -A D T lp C p -B T 2 L p )Xcl\} ,(107) 

= ^{(d\ 2 D x 2 G^C 2 ~b\q^J})Lc\+{d\ 2 D 12 G C 2 -D* 2 C l -B 2 Qoo)L^C 2 +(AD\ p D lp G^C 2 

-B T 2 Lf)XC T 2 +a D T ip D lp GC 2 -A D T lp C p -BlL p )X^cl] . (108) 
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Expressions for L ^ and obtained by differentiating equations (57) and (58), are given by: 

0 = A T L^LfA* [A« %+Z, A lJ \MC T p C p ) (fl ] . (109) 

O^Ax^+^A^lA^X+xA ( i ,T + (B p S T p ) fil ] . (110) 

while L (J) , Q ( J\ and j?°Ve the same as the H „ case, equations (100), (101), and (77), respectively. 
From equation (52), 

(C;e,)« = (-5 lp o0 1 C 2 ) r (C,-B„CC 2 )+(C f -D, p GC 2 ) :r (-5 lp G w C 2 ) . (Ill) 

2. Computation of H a . As with the computation of the hessian, the derivative of the 
homotopy function with respect to a is given by: 

= — {2f(DpD 12 GC r Df 2 C l -5 2 r a.)l-C 2 +(A5f p 5 1() GC 2 -ADf J) C p -S 2 L p )X c[]} ,(1 12) 
da da 

= 2[{t\ 2 D \ 2 GC 2 +D 12 GC 2 -t T l2 C \-B 2 (f*)LC 2 +(D J 2 D \ 2 GC 2 

-D[ 2 C 1 -fi[(2J^C 2 ^(AD[^ lp GC 2 -AD^C p -BlVXc[] . (113) 

Expressions for L p and X are obtained by differentiating equations (57) and (58) to obtain 

• (ll4) 

0 = Ax+xA t hax+xA t +b i ,bI+B i ,b t p ] , (115) 

and 

^ = ^max~Anin • (116) 

Note that the procedure described in section III. A, where the H „ and H 2 homotopies are 
performed distinctly, simplifies the computations significantly in that the plant matrices remain fixed 
and only yor A are varied at one time. 


IV. DESIGN EXAMPLE 


To demonstrate the homotopy algorithm applied to optimal controller synthesis, the four-disk 
example originally described in reference 13 and more recently by numerous others (see reference 
11) will be used. The four-disk model used in the example problem was derived from a laboratory 
experiment and represents an apparatus developed for testing of pointing control systems for flexible 
space structures with noneolocatcd sensors and actuators. As illustrated in figure 5, tour disks are 
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Ml, 01 


M2, 02 


M3, 03 


M4.04 



12 3 4 

Figure 5. Four-disk system. 


rigidly attached to a flexible axial shaft with control torque applied to selected disks and the angular 
displacement of selected disks measured. The equations of motion may be written as: 


/j0 !+*!«?! -0 2 ) = O , 

(117) 

/ 2 0 2+^2 ( 0 2 — 0 3 ) — ^1 ( $1 ~@2 ) = 0 , 

(118) 

/ 3 0 3 +^ 3 ( 6)3 —$4 )—K 2 ( $2 — 03 ) = 0 , 

(119) 

i 4 e 4 -K 3 (d 3 -e 4 ) = o , 

( 120 ) 

or 


Mc}+Kq = Bu , 

( 121 ) 

where the generalized displacements are the angular displacements of the disks, q 
and the input vector consists of the moments applied to each disk, u T = [M x M 2 M 3 
state vector as x T = [q T q 7 ] results in the state space formulation: 

7 =[01 02 03 0 4 ]> 

M 4 ], Defining the 

x = Ax + Bu , 

( 122 ) 


where 


A = 

0 I 

B = 

0 


M~ X K 0 


M~ ] R 


(123) 


For simplicity, the stiffness and inertia terms are set to unity ((G7/L), = K ( = /, = 1 , i=l:4). In this 
case, the mass matrix is a 4 by 4 identity matrix and the stiffness matrix is: 


K = K-K = K 


1-10 0 
-1 2-1 0 
0-1 2-1 
0 0-11 


(124) 
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The plant is modeled with parametric uncertainty corresponding to uncertainty in the stiffness of each 
shaft. Let uncertainty in the shaft stiffness be modeled as: 

K = Kq+AK => K= (Kq+AK)R . (125) 


Hence, the A matrix becomes: 


A = 


0 I 
£ 1-1 


+ 


o 

= Aq+AA . 


0 

-m _1 £ 


AK[I 0] , 


A block diagram representation of the plant with uncertain stiffness is shown in figure 6. 



(126) 

(127) 


Figure 6. Plant with uncertain stiffness. 


A. The H 2 Case 

To demonstrate the homotopy algorithm applied to H 2 controller synthesis, the same four- 
disk model used in reference 11 will be used here. These results provide a direct comparison 
between the homotopy algorithm of reference 1 1 (HAO) and the algorithm presented in section III 
(H2HOM). The main distinction between the two homotopy algorithms for the H 2 case is the 
compensator architecture. HAO employs a general architecture that may be restricted to various 
parameterizations including the controllability canonical form (the controller canonical form is used in 
H2HOM). When the controllability canonical form is implemented in HAO, as in this example, the 
compensator is still represented in a general architecture, resulting in the evaluation of five 
necessary conditions (equation (21)). The HAO code has been highly optimized for efficient 
computation with the result that superfluous computations are not evaluated. The homotopy 
algorithm of this paper, H2HOM, is patterned after the general approach of HAO and utilizes some 
of the more efficient computational aspects of the HAO code. 

The control design philosophy for this example is to scale the nominal control weight and the 
nominal sensor noise intensity by the parameter q. As q is reduced, the control authority is 
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increased. For comparison with the result published in reference 11, a full-order (eighth-order) 
compensator was synthesized. Although the results can be directly obtained from the LQG Riccati 
equations, the full-order compensator was chosen to tax the H2HOM algorithm, which must 
optimize over a greater number of parameters with increasing compensator order. 

Table 1 shows a comparison of the results from the H2HOM and HAO algorithms for the full- 
order compensator along with results from the H2HOM algorithm for sixth- and second-order 
compensators. All pertinent parameters as well as logic for step size scaling and the computation of 
the prediction and correction errors are identical in both algorithms, which are implemented in 
MATLAB™ on a 486 66MHz computer. Whereas the HAO code required a minimum step size of 
1.907e-7, the H2HOM code was much better conditioned and required a minimum step size of 0.025. 
As a consequence of the smaller step sizes with HAO, 2,504 hessian computations were required as 
opposed to only 63 hessian computations with H2HOM. The HAO code has been tuned extensively 
for efficient computation as is reflected in the small number of flops required. In spite of the 
significantly smaller number of flops required for HAO, the H2HOM code required significantly less 
clock time for convergence to the same final compensator. (The results generated by the author using 
the HAO code differ slightly from those reported in reference 11, although the parameters in the 
HAO algorithm are the same. It is likely that the published results were generated with an earlier 
version of the HAO code. The qualitative trends remain the same.) 

Table 1. Comparison of H2HOM and HAO algorithms. 


Algorithm 

HAO 

H2HOM 

H2HOM 

H2HOM 

Compensator Order 

8 

8 

6 

2 

Number Hessian Comp. 

2,504 

63 

60 

30 

Minimum Step Size 

1.907e-7 

0.025 

0.05 

0.1 

Maximum Step Size 

0.1 

0.1 

0.1 

0.2 

Max. Number Correction Iterations 

9 

7 

8 

5 

Mflops 

287 

936 

455 

37 

Time (s) 

5,104 

883 

488 

73 


In reference 11, the controllability canonical form is assessed as poorly conditioned because 
of the small minimum step size. However, table 1 indicates that the static gain formulation in 
H2HOM yields a substantial improvement in conditioning along the homotopy path over the HAO 
implementation of the canonical compensator. However, this may not be the case in general due to 
the tendency toward ill-conditioning characteristic of canonical forms. An even more significant 
benefit of this formulation is the straightforward extension to the and mixed problems. 


B. The Mixed HilH^ Case 

The homotopy algorithm was also used to synthesize reduced-order mixed 
compensators for the four-disk problem. The problem was formulated in terms of a robustness 
design and a performance design. Similar to the methodology of Luke, 14 the portion uses 
weighted sensitivity for the tracking problem and minimizes the control energy due to stochastic 
disturbances using the // 2 norm. The block diagram for this problem is shown in figure 7. The output 
is angular position of disk 3, and the control input is torque applied to disk 3. 
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Ze dml dm2 Zu da d 



Figure 7. Closed-loop four-disk block diagram. 

First, a nominal performance design was performed by varying the sensitivity weight. We, to 
obtain desirable step responses. Then as described in section IV, the additive uncertainty 
corresponding to uncertain shaft stiffness was included. This uncertainty representation is somewhat 
conservative due to the complex uncertainty representation, but still serves to allow robustness to 
uncertainty in the modes, which in turn adds to the stability margins. The uncertainty weights were 
then scaled to allow the infinity norm of the closed-loop system to be less than 1. For the //„ 
minimization problem, the input vector w consisted of four inputs corresponding to the uncertainty 
model (fig. 6); measurement noise, dm\ \ disk 3 position command, 03 ,com'> and actuator noise, da. The 
output vector z consisted of four outputs corresponding to the uncertainty model; control energy, Zw; 
and weighted error, Ze. For the H 2 minimization problem, the input vector w p consisted of a 
stochastic disturbance torque, d, and measurement noise, dm2. The output p was the control energy 
Z u . The sensitivity weight was given by 


A full-order compensator was synthesized for the portion using the standard 2-Riccati 
(DGKF) solution 15 for comparison with the full- and reduced-order compensators synthesized by 
the homotopy algorithm. The homotopy algorithm was able to reproduce the 2-Riccati solution with 
the exception that the code was stopped prior to obtaining the minimum gamma since convergence 
slows considerably in the neighborhood of optimum. For the sixth-order compensator, the 
homotopy was terminated at y = 1.08. Figure 8 shows the maximum singular value plots of the 
closed loop for each H «, compensator. It should be emphasized that using schur-balanced model 
reduction and optimal hankel norm model reduction on the full-order H compensator, it was not 
possible to obtain a stabilizing reduced-order compensator of any order. Thus, the homotopy 
algorithm was able to synthesize fixed-order H <*, controllers where standard model reduction 
techniques failed. 

H„ compensators arc generated as a special case of the mixed H 2 IH „ homotopy with the H 2 
norm weight A set to zero. Then for fixed y, a homotopy is performed on A, which is increased and 
the minimum H 2 norm is obtained. Simultaneously, the norm increases and the actual infinity 
norm approaches the overbound. This is the key benefit of the mixed H 2 IH « formulation: the 
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Figure 8. Singular values for H«, compensators. 


conservatism of the design is reduced by reducing the gap between the overbound y and the 
actual Ho* norm, while the performance increases by minimizing the H 2 norm. 


The improvement in performance as indicated by the lower H 2 norm is shown in figure 9, 
where the minimum H 2 norm is plotted versus the //«, norm as the overbound y is decreased. For the 
DGKF full-order solution, note that the largest //« norm attainable is 1.7, which limits the attainable 
H 2 norm to 1.088. However, using the mixed H 2 /Hoo formulation, the homotopy algorithm generates a 
full-order compensator for a wider range of H„ norms over which smaller H 2 norms are attainable for 
a given y overbound. In this case, the largest //„ norm was 3.795 with an H 2 norm of 0.3317, which 
is a 70-percent reduction of the solution. Similarly with the sixth-order mixed H 2 !Hoo 
compensator, the minimum H 2 norm attained was 0.52 with an Hoc norm of 4.93. Thus, the presence 
of the gap between the overbound and the actual norm limits the performance range attainable 
with the Hoo solution, but is removed when the mixed H 2 !Hoo formulation is employed. 



H— Infinity Norm 


Figure 9. H „ versus H 2 norm trades. 
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V. CONCLUSIONS 


A novel homotopy algorithm is developed to synthesize fixed-order H 2 and compensators 
employing a controller canonical form, and a representative flexible structure is used to demonstrate 
the numerical results. These results indicate that the static gain optimization formulation may be a 
more efficient means of synthesizing dynamic compensators than employing a general compensator 
architecture. An even more significant benefit of this approach is the straightforward extension to the 
fixed-order //„ and mixed H^IH^ control problems. The synthesized reduced-order compensators 
perform well when compared to full-order controllers, which is highlighted by the fact that standard 
controller order reduction techniques do not yield a stabilizing compensator. The fixed-order mixed 
H2IH00 formulation is shown to offer improved performance over standard H „ compensators by 
minimizing the H 2 norm while removing (or reducing) the gap between the actual H„ norm and the 
gamma overbound. 
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